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Weak lensing of the CMB changes the unlensed temperature anisotropy and polarization power 
spectra. Accounting for the lensing effect will be crucial to obtain accurate parameter constraints 
from sensitive CMB observations. Methods for computing the lensed power spectra using a low- 
order perturbative expansion are not good enough for percent-level accuracy. Non-perturbative 
flat-sky methods are more accurate, but curvature effects change the spectra at the 0.3-1% level. 
We describe a new, accurate and fast, full-sky correlation- function method for computing the lensing 
effect on CMB power spectra to better than 0.1% at I < 2500 (within the approximation that the 
" " ", lensing potential is linear and Gaussian). We also discuss the effect of non-linear evolution of the 

^— ^ ■ gravitational potential on the lensed power spectra. Our fast numerical code is publicly available. 
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I. INTRODUCTION 

The CMB temperature and polarization anisotropies are being measured with ever increasing precision. The 
statistics of the anisotropies already provide valuable limits on cosmological parameters, as well as constraints on 
early-universe physics. As we enter the era of precision measurement, with signal-dominated observations out to 
small angular scales, non-linear effects will become increasingly important. One of the most significant of these over 
scales of most interest for parameter estimation is weak gravitational lensing by large scale structure. Fortunately it 
can be modelled accurately as a second-order effect: the linear gravitational potential along the line of sight lenses the 
linear perturbations at the last scattering surface (see e.g. Refs. 00,0 and references therein). Modelling of fully 
, non-linear evolution is not required for the near future on scales of several arcminutes (corresponding to multipoles 
I < 2000) for the temperature and electric polarization power spectra. Non-linear corrections can easily be applied to 
] the lensing potential if and when required, provided that its non-Gaussianity can be ignored 

In principle, the weak-lensing contribution to the observed sky can probably be subtracted given sufficiently accurate 
and clean high-resolution observations. Early work in this area la, El, suggested a limit on the accuracy of this 
reconstruction due to the statistical nature of the (unknown) unlensed CMB fields. More recently, it has been argued 
that polarization removes this limit in models where lensing is the only source of i?-mode polarization on small 
scales . If subtraction could be done exactly we could recover the unlensed Gaussian sky, and use this for all further 
t-^ ■ analysis. However current methods for subtracting the lensing contribution are approximate, and not easy to apply 
to realistic survey geometries. The result of imperfect lensing subtraction is a sky with complicated, non-Gaussian 
statistics of the signal, and significantly more complicated noise properties than the original (lensed) observations. 
For observations in the near future, a much simpler method to account for the lensing effect is to work with the 
lensed sky itself, modelling the lensing effect by the expected change in the power spectra and their covariances. The 
effects of lensing non-Gaussianities on the covariance of the temperature and i?-mode polarization power spectra are 
likely to be small, but this will not be the case for the J3-mode spectrum once thermal-noise levels permit imaging 
of the lens-induced B modes 9j. In this paper we discuss how to compute the lensed power spectra accurately. The 
simulation of lensed skies and the effect on parameter estimation is discussed in Ref. |l(j • 

On scales where the non-Gaussianity of the lensing potential can be ignored, the calculation of the lensed power 
spectra is straightforward in principle. However, achieving good accuracy on both large and small scales for all the 
CMB observables is surprisingly difficult. The lensing action on the CMB fields at scales approaching the r.m.s. of 
the lensing deflection angle (~ 3 arcmin) cannot be accurately described with a first-order Taylor expansion, as in the 
full-sky harmonic method of Ref. 0. There is not much power in the unlensed CMB on such scales, but a first-order 
Taylor expansion still gives lensed power spectra that are inaccurate at the percent level for I > 1000. The lensed 
CMB on scales well below the diffusion scale is generated by the action of small-scale weak lenses on the (relatively) 
large-scale unlensed CMB, and a Taylor expansion should become more accurate again 11]. (However, non-linear 
effects are also important on such scales.) The breakdown of the Taylor expansion can be easily fixed by using the 
flat-sky correlation- function methods of Refs. 00, which can handle the dominant effect of the lensing displacement 
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in a non-perturbative manner. However, a new problem then arises on scales where the flat-sky approximation is not 
valid. As noted in Ref. Q , this is not confined to large scales due to the mode-coupling nature of lensing: degree-scale 
lenses contribute significantly to the lensed power over a wide range of observed scales. In this paper we develop 
a new method for computing the lensed power spectra that is accurate on all scales where non-Gaussianity due to 
non-linear effects is not important. We do this by calculating the lensed correlation functions on the spherical sky. 
This allows us to include both the non-perturbative effects of displacing small-scale CMB fluctuations, and the effects 
of sky curvature. 

This paper is arranged as follows. We start in Section[H]with a brief introduction to CMB lensing, then in Section lTTTl 
we review previous work on flat-sky correlation-function methods and present our new full-sky method and results. In 
Section lTVl we compare our new results with the flat-sky correlation- function results of Refs. [HQ and the perturbative 
harmonic result of Ref. and explain why the latter is not accurate enough for precision cosmology. The effect 
of non-linear evolution of the density field on the lensed power spectra is considered in Section [V] We end with 
conclusions, and include some technical results in the appendices. 



II. CMB LENSING 



Gradients in the gravitational potential transverse to the line of sight to the last scattering surface cause deviations 
in the photon propagation, so that points in a direction n actually come from points on the last scattering surface in 
a displaced direction n'. Denote the lensed CMB temperature by 9(n) and the unlensed temperature by O(n), so 
the lensed field is given by 0(n) = O(n'). The change in direction on the sky can be described by a displacement 
vector field a(n) = V-0, so that (symbolically) n' = n + Vi/>. Here ip is the lensing potential which encapsulates the 
deviations caused by potentials along the line of sight. More rigorously, on a unit sphere the point n' is obtained from 
n by moving a distance |V-0| along a geodesic in the direction of Vi/'(n), where V is the covariant derivative on the 
sphere |12( . We assume that the lensing is weak, so that the potentials may be evaluated along the unperturbed path 
(i.e. we use the Born approximation). Lensing deflections are a few arcminutes, but are coherent over degree scales, 
so this is a good approximation. 

In terms of the zero-shear acceleration potential ty, the lensing potential in a flat universe with recombination at 
conformal distance x* is given by the line-of-sight integral 

t/)(n) = -2 dx*(xn;77o-x)^^. (1) 
Jo XX* 

Here we neglect the very small effect of late-time sources, including reionization, and approximate recombination as 
instantaneous so that the CMB is described by a single source plane at x — X* ■ The quantity r/o — \ is the conformal 
time at which the photon was at position %n. With the Fourier convention 

*( x ^)=/(^*( k ^) e * X ' ( 2 ) 



and power spectrum 



(*(k; ry)**(k'; rf)) = |^P*(fc; rj, r/)8(k - k'), (3) 



the angular power spectrum of the lensing potential ip evaluates to 

Cf =167r / f F'dx [ X *a^V*(k; m -x,Va-X%(kx)3i(kx')(^pt) (^f) - (4) 
J K Jo Jo V X*X J \ X*X J 

In linear theory we can define a transfer function T^,(k,rj) so that ^(k;^) = T^{k, rj)lZ(k) where 7Z(k) is the 
primordial comoving curvature perturbation (or other variable for isocurvature modes). We then have 



Cf = Wir J ^ V n {k) £' d X T*(k; Vo - x)ji(kx) (^jr/) 



(5) 



where the primordial power spectrum is Vn(k). This can be computed easily numerically using CAMB 1 |13|. and a 
typical spectrum is shown in Fig. ^ 
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FIG. 1: The power spectrum of the lensing potential for a concordance ACDM model. The linear theory spectrum (solid) is 
compared with the same model including non-linear corrections (dashed) from halofit pH using Eq. 



III. LENSED CORRELATION FUNCTION 



A. Flat-sky limit 

We start by calculating the lensed correlation function in the flat-sky limit, broadly following the method of Ref. pj. 
We use a 2D Fourier transform of the temperature field 



e(x) = J fjeay 1 -*, 



and the power spectrum for a statistically isotropic field is then 

(e(i)e*(i')> = cp5(i-i'). 

Lensing re-maps the temperature according to 

6(x) = 6(x + a), 



(6) 



(7) 



(8) 



where in linear theory the displacement vector a is a Gaussian field. We shall require its correlation tensor 
(aj(x)oj(x')) to compute the lensed CMB power spectrum. Introducing the Fourier transform of the lensing po- 
tential, ip(l), we have 



a(x)=i J ^lV-(l)e ll x , 



so that 



(ai(x)aj(x')) = 



d 2 l 



(27T) 



LLC?e a < x - x ">. 



(9) 



(10) 



By symmetry, the correlator can only depend on Sij and the trace- free tensor r/j7j\, where r = x — x'. Evaluating 
the coefficients of these two terms by taking the trace of the correlator, and its contraction with rV J , we find 



(o,-(x)ai(x')) - i- f dll 3 Cfjo(lr)5 i:j - i- / dl l 3 Cf J 2 (h !/ , , , . 



(11) 
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where J n (x) is a Bessel function of order n. Note that the trace-free term is analytic at r = due to the small- r 
behaviour of J^ilr). Following Ref. 1J, let us denote (a(x) • a(x')) by C g i(r) so that 



C g lW = ^ I nllUr-y (12) 



Similarly we define the anisotropic coefficient 



so that 



C gl , 2 (r) = i- / dZZ 3 Cf J 2 (/r). ,1.:-!) 



(ai(x)aj(x')) = ^C g i(r)% - C g i, 2 (r)f {i f j} . (14) 



The lensed correlation function £(r) is given by 



£(r) = (6(x)6(x')) 

-Cfe il - r (e il '[ a; W- Q! ( x '^}, (15) 



d 2 l 

{2nf 



where we have assumed that the CMB and lensing potential are independent (i.e. we are neglecting the large scale 
correlation that arises from the integrated-Sachs- Wolfe effect and has only a tiny effect on the lensed CMB). Since 
we are assuming a is a Gaussian field, 1 • [c*(x) — ct(x')] is a Gaussian variate and the expectation value in Eq. 115|) 
reduces to 

^l-Kx)-*^')]) = exp ^— i<[l • (a -ex')} 2 ) 

= exp (^-izV(r) +cos2(0! - </> r )C gl>2 (r)]) , (16) 

where we have used rPf^fj) — Z 2 cos2(</>i — </> r )/2 and defined <r 2 (r) = C g i(0) — C g \(r). Here, e.g. (f>\ is the angle 
between 1 and the a;-axis. The cos2(0i — <j) r ) term in Eq. (|16fl is difficult to handle analytically. Instead, we expand 
the exponential and integrate term by term. Expanding to second order in C g i i2 , we find 



* (r) = h / w ^ e ~' 2CT2(r)/2 [( 1 + iV c ^ 



(17) 



Expanding to this order is sufficient to get the lensed power spectrum to second order in Cf; higher order terms in 
C g i, 2 only contribute at the O(10~ 4 ) level on the scales of interest. Note that the exp(— l 2 a 2 /2) term is easily handled 
without resorting to a perturbative expansion in Cf . Since a 2 is significantly less than C g i i2 (as shown in Fig.[2J), the 
perturbative expansion in C g i i2 converges much faster than one in a 2 . Equation l|17f) extends the result of Ref. [I| to 
second order in C g i 2 . 



1. Polarization 

The polarization calculation is also straightforward in the flat-sky limit llj. We use the spin —2 polarization 
P = Q + ill, where Q and U are the Stokes' parameters measured with respect to the fixed basis composed of the x 
and — y axes. Expanding P(x) in terms of the Fourier transforms of its electric (P) and magnetic (B) parts, we have 

P(X) = " / S {E{1) ~ ^(OK^e* 1 *, ( 18 ) 

where (d x — id y ) 2 e lhx /l 2 = — e - 2i </ , i e llx [ s a S pi n _2 flat-sky harmonic. The polarization correlation functions are 
defined as 

e+(r) = (e^'P*(x)e 2 *P(x')), (19) 
£_(r) = (e 2i ^P(x)e 2i< ^P(x')}, (20) 
= <e(x)e 2 *P(x')), (21) 
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FIG. 2: The functions cr 2 (/3) = C g i(0) - C g i(/3) [solid] and Cgi, 2 (/3) [dashed] as a function of angular separation /3 (in radians) 
for a typical concordance model. The results are calculated using the full-sky definitions of Eqs. (13511 . and use the linear power 
spectrum for Cf . 



where TT — cj> r is the angle to rotate the x-axis onto the vector joining x and x', so that e.g. e 2l ^ r P(x) is the polarization 
on the basis adapted to x and x'. Then the lensed correlation functions to second-order in C g i t2 are 



1 

2^ 



1 

27 



ZdZ (Cf + Cf)e- iV2 W/2 



Zdi (Cf - Cf )e-' V2 W/ 2 



16 



* 4 C£» Jo(Zr) 



+ \l 2 C g , a {r)J 2 {lr) + ±-l^ 2 (r).h(lr) 



1 



l + y/C g \ 2 (r)j J 4 (Zr) 
+ \l 2 C^ 2 {r) l -[J 2 {lr) + J 6 (lr)] + ±-l*C 2 glt2 (r)±[J (lr) + J 8 (lr)} 
1 + -^ 4 Cgl>)) J 2 (Zr) + I/ 2 C g , 2 (r)i[J (Zr) + J 4 (Zr)] 
+ T^^ g 2 i,2W^[J2(/r) + J 6 ar)] 



(22) 



-(23) 



(24) 



Here Cf and Cf are the -E-mode and £?-mode power spectra, and Cf is the Q-E cross-correlation. This is the 
straightforward extension of the result in Ref. 01 to higher order; 2 see that paper for further details of the calculation. 

The lensed £+(r) has the same structure as for the temperature since the unlensed correlation functions involve 
the same Jo (Ir) , and there are no complications due to the different local bases defined by the displacement r and its 
image under lensing r — a' + a since the phase factors from the rotations cancel. This is not the case for the lensed 
|_(r) and £x(r). 



Note that we disagree with the statement in Ref. 1 141 that a 0(cf) expansion is very accurate. Indeed CMBFAST 4.5 actually uses the 
non-perturbative <r 2 term (as advocated here) rather than the lowest-order series expansion given in the paper. 
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2. Limber approximation 

At high I the power spectrum V^(k) varies slowly compared to the spherical Bessel functions in Eq. which pick 
out the scale k ~ l/ X - Using 

" k^kMkxMkx') = ^S( X ~ X% (25) 



2x 2 



we can Limber-approximate Cf as 



c t^^~ I*' xd x v*(i/x;m-x)(^^) 2 ■ (26) 



i d Jo V x*x 

Changing variables to fc = i/x, we find 



C gl (r) w 4^ y dfc y d X P*(fc;77o - x) (^T^) J o( fc Xr), (27) 

in agreement with Ref. (ij if we note that his P$(k) — 7 ? *(A:)/(47rfc 3 ) outside radiation-domination. Ref. also 
defines C g \^{r) as (in our notation) 

C gl , 2 (r) ^Anjdk J dxr*(k;r, - X ) (~ -M^r), (28) 

which is the Limber-approximation version of Eq. For the results of this paper we do not use the Limber 

approximation, though the approximation is rather good. 



B. Spherical sky 

The fiat-sky result for the lensed correlation function is non-perturbative in <J 2 (r) and this turns out to be crucial 
for getting high accuracy in the lensed power spectrum on arcminute scales. Consider the contribution to the lensed 
correlation functions from the unlensed CMB at multipole I. Both a 2 and C g i 2 appear with a factor I 2 and so 
the (dominant) I 2 a 2 term cannot be handled accurately with a low-order expansion at high I. Physically, this is 
because the typical lensing displacement is then comparable to the wavelength of the unlensed fluctuation, and so 
approximating the fluctuation as a gradient over the scale of the lensing displacement is inaccurate. The error from 
this gradient approximation on the lensed power spectra will be large on any scale |1| where the dominant contribution 
is from unlensed fluctuations with wavenumber 1' comparable to the typical lensing displacement at scale |1 — l'|. 

As noted in the Introduction, the small-scale cut-off in the power in the unlensed fluctuations due to diffusion 
damping means that the gradient approximation should not get uniformly worse on small scales; the approximation 
should be poorest on scales of a few arcminutes. We also noted that the the flat-sky approximation will be suspect on 
large scales, and also on any scale where the dominant contribution is from large-scale lenses, i.e. those for which their 
mode-coupled wavenumber |1 — 1'| small. What is needed for an accurate calculation (on all scales where non-linearities 
in the lensing potential are not important), is a non-perturbative treatment of a 2 and a proper treatment of curvature 
effects in the correlation functions. In this section we show how to generalize the flat-sky calculation to spherical 
correlation functions. 

On the full sky we can expand the temperature field in spherical harmonics 

9(n) = ^9 ;m ^ m (n), (29) 

Irn 

and the temperature correlation function is defined by 

m = (9(^)6(112)), (30) 

where (3 is the angle between the two directions (rii • A2 = cos 0) . The power spectrum is defined as the variance of 
the harmonic coefficients Cf = (|Oz m | 2 ) for a statistically-isotropic ensemble. 
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FIG. 3: The geometry of the weak lensing deflections (shown without curvature for clarity). 



We define a spin-1 deflection field \u = a ■ (eg + ie^), where eg and are the unit basis vectors of a spherical-polar 
coordinate system. Rotating to the basis defined by the geodesic connecting rii and n2, the spin-1 deflection (denoted 
with an overbar in the geodesic basis) has real and imaginary components 

ai cosi/h = 5?ia(ni), ot\ sin tpx — 3ia(ni), (31) 

and similarly at h 2 . Here, ol\ — |a(ni)| is the length of the lensing displacement at fii and ipi is the angle it makes 
with the geodesic from fix to h 2 (see Fig- EJ ■ In terms of these angles we have the lensed correlation function 

£(/3) = <e(ni)6(n' 2 )> (32) 
= (33) 

Int 

= 52 Cf4m'C9)Wt n («i,^i)^/(a2,^)}. (34) 

Irara' 

The easiest way to see the last step is to put ni along the z-axis, and \\i in the x-z plane so that has polar coordinates 
(ai, ipi). The harmonic at the deflected position h' 2 can be evaluated by rotation: YJ m (n 2 ) = [l) _1 (0, (3, 0)F; m ](a2, 1P2), 
where [l)Y; m ](n) is a spherical harmonic rotated by the indicated Eulcr angles. We have neglected the small correlation 
between the deflection angle and the temperature so that they may be treated as independent fields. The remaining 
average is over possible realizations of the lensing held. 

We assume the lensing potential is Gaussian, so the covariance of the spin-1 deflection field can be determined using 
the results 

0/ _L 1 

(Mni)Mh 2 )) = -52-£-Ki + i)c+dL u {p) = -CM(P), 

I 

2/ 1 I 

( 1 a*(n 1 )ia(n 2 )} = ]T + l)Cf d^fi) = C gl (/3). (35) 

I 

As in the flat-sky limit, it is convenient to define er 2 (/3) = C g i(0) — Cgi(j3). The covariance of the Gaussian variates 
3?icv(n 1 ), 3 1 a(n 1 ), Kia(n 2 ) and 3ia(n 2 ) are determined by Eq. (|35|l . Transforming variables to a%, tpi, a 2 and ip 2 
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we find their probability distribution function 

4aia 2 e 2 

Pr(ai,a2, ipi,ip2) 



(ai cos ipi+a 2 cos j/) 2 ) / (t +2C g i— C g i, 2 ) 



(2tt) 2 + 2C gi " C gl,2 

— 5(0:1 sini/'i+c.2 sini/j 2 ) 2 /(o' 2 +2C g i+o g i i 2; — 2 



1 (a! sini/'i+C!2 sini/i 2 ) 2 /(o" 2 +2C gl +C g i : 2) — cos^i-q 2 cos i/j 2 ) 2 /(o ,2 +C g i, 2 ) 



y/<7 2 + 2C* g l + C g l, 2 V^+CgU 
p -2"("i sinVi-C( 2 sin V , 2) 2 /(c 2 -C' gli2 ) 

x 7t— ■ ( 36 ) 

Here and below we have left the dependence of a 2 , C g i and Cgj 2 on /3 implicit. Our general strategy to evaluate 
Eq. (|34[1 is to expand Pr(ai, 012, tpi, 1P2) in Cgi and Cgi^, but not a 2 , before performing the integral over the angles ipi 
and tp2 in the expectation value. The remaining integrals over a% and 012 then enter through functions of the form 

r°° 2rr / n \ i 2,2 
X imn ^J Q ^ (^) e~ a d l mn (a)da. (37) 

Since terms involving C g i(/3) are suppressed at high I (they do not appear in the flat-sky results), while at low I the 
leading-order result neglecting C g i and C g i !2 altogether is very accurate, we neglect terms involving C g i entirely. This 
approximation is very accurate [< O(10~ 4 )] (for completeness the full second-order result is given in Appendix [UJ. 
As shown in Fig. the values of C g i !2 are much smaller than tr 2 , so a perturbative treatment in C g \ t 2 is sufficient as 
in the flat-sky case. Working to second-order in C s i t 2, we find 

£ ~ £ 4^ Q e | ^"000^00 + — ^ Cgi, 2-^000^1-1 + Cgi, 2 (^ockWoo + -^220^2-2) j> (38) 

where primes denote differentiation with respect to a 2 [note that the Xi mn are implicit functions of /3 via the de- 
pendence on cr 2 (/3)]. In Appendix iBl we develop approximations for the integrals Xi mn which are accurate for all I. 
Applying these approximations, the required Xi mn are 

Aooo « e- l ^° 2 l* (39) 
A 220 w \ V(l + m - mi + l) e - [i( ' +1) - 2lff2/4 . (40) 

The expansion of these results to 0(a 2 ) may also be derived straightforwardly by using the series expansion of d\ nn {a) 
for small a. (The smallness of a 2 guarantees that the integral is dominated by the small a region). However, it is 
important to retain the correct non-perturbative form for high I. 

In the limit of large I the limiting result d l mn (/3) — > (— I)"~ m J m _„(//3) shows that the full result of Eq. 1(55)1 reduces 
to Eq. I|17|) in the flat-sky limit and is therefore consistent. In the limit in which the separation angle (3 — > we have 

1 1 

where C? is the lensed power spectrum. This expresses the fact that weak lcnsing does not change the total fluctuation 
power. 



1. Polarization 



We can extend the previous calculation to polarization. Defining Stokes' parameters with the local x-axis along the 
^-direction and y along —<p, the quantities Q ± iU are spin =p2 respectively. We can expand Q ± iU in terms of the 
spin-weight harmonics as [15( 

(Q ± iU)(n) = J2( E im T iBi m ) T2 Yi m (n), (42) 

irn 

which expresses P = Q + iU as the sum of its electric (E or gradient-like) and magnetic (B or curl-like) parts. (Our 
conventions for the polarization harmonics and correlation functions follow Refs. [To . Il7| ; see these papers for a more 
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thorough introduction). The polarization correlation functions can be defined in terms of the spin ±2 polarization 
defined in the physical basis of the geodesic connecting the two positions. As for the temperature, we evaluate the 
polarization correlation functions by taking rii along the z-axis and n 2 in the x-z plane at angle j3 to the z-axis. 
With this geometry, the polar-coordinate basis is already the geodesic basis connecting fix and n 2 so that the lensed 
correlation functions are 

£+(/?) = (P*(fii)P(fi a )), (43) 
£_(/3) = (P(n 1 )P(n 2 )}, (44) 
lx{(3) ee (e(n 1 )P(n 2 )). (45) 

Under a lensing deflection the polarization orientation is preserved relative to the direction of the deflection (we are 
neglecting the small effect of field rotation 0), i.e. the polarization undergoes parallel transport. The geometry of 
the deflections is shown in Fig. [21 

We can easily evaluate the lensed polarization on the connecting geodesic basis (between fii and n 2 ) as 

P(ni) = Picture- 2 ** 1 . (46) 

The rotation angle is that needed to rotate the spin —2 polarization from polar coordinates (coinciding with the 
the ni-n'j basis at n^) to the geodesic basis connecting ni and n 2 . 

For the lensed polarization at n 2 a little more work is required. Let x' denote the angle between the geodesies 
connecting n 2 to n' 2 , and ni (along the z-axis) to n' 2 (see Fig.|2Jl. The lensed polarization at n 2 on the geodesic basis 
adapted to ni and n 2 is then 

P(n 2 ) = P(ii' 2 )e 2lx ' e- 2 ^ 2 . (47) 

We can write n' 2 as the direction obtained by rotating a direction with polar angles (a 2 , 1/^) by an angle (3 about the 
y-axis, i.e. n 2 = D(0, /3, 0)(a 2 , ^2)- Writing P as (Q — ill)*, and using Eq. @2J, we have 

P(h) = ^(E lm + »B Im )* +a y^(n). (48) 

I in 

Using the rotation properties of the spin-s harmonics (see Appendix [XJ , we then find 

P(n 2 ) = e ^x' e -2i^ e -2i K J- (£ im +iB| m )*DL'(0,A0)2C(«2,^). (49) 

Imm' 

The angle n is the rotation about n 2 that is required to bring the polar basis there onto that obtained by rotating 
the polar basis at (a 2 ,'0 2 ) with D(Q,(3, 0). Since the latter is aligned with the geodesic basis adapted to n 2 and n' 2 , 
we have K = x' an d the lensed polarization at n 2 simplifies to 

P(n 2 ) = e" 2 * Yl (Eim+ l B lm rd l mm ,((3) 2 Y l * m ,(a 2 ^ 2 ). (50) 

Imm' 

We can now quickly proceed to the following expressions for the lensed polarization correlation functions: 

f+(/3) = £(Cf + Cf)4 im '(/3)(e 2 *2^ m (ai,^i)2^(a2,^)e- 2 ^), (51) 

Imm' 

= ^(Cf -Cf)4 im ,(/3)(e- 2 ^_ 2 ^ m (a 1 ,^ 1 ) 2 y ; * m ,( a2 ,^2) e - 2 ^), (52) 

Imm' 

ixW) = E Q X rfU(/3)Wm(« 1 ^ 1 ) 2 F4'(«2,V 2 )e- 2 *), (53) 

/mm' 

where the expectation values are over lensing realizations. Here, and C; B are the power spectra (\Ei m \ 2 ) and 
(|-Bzm| 2 ) respectively. The cross-correlation power spectrum is C ; x = (Qi m E* m ). 

We evaluate the expectation values in Eqs. (|51H53() following the earlier calculation for the temperature, i.e. ex- 
panding Pr(ai, a 2 , "01, ^2) to second order in C g i j2 before integrating. As for the temperature, C g i terms contribute 
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negligibly (see Appendix IU1 for the full result). We find the following results for the lensed polarization correlation 
functions to second order in C g \ >2 '- 



£+ w 4?r (Cf + C^) {-^022^22 + 2(^1,2X132X121 d 31 + Cg l2 [(X' 022 ) 2 d l 22 + X 2i2 X 220 d l 40 ]^ , (54) 

Imm 



I- ~ X/ ~TZ~(Cf ~ Cf)\ Xl 22 d l 2 _ 2 + C g l,2[X!2idi _i + X! 32 c?3_3 



Air 

Imm' 



2 < -'gi,2[2(XQ 22 ) 2 (i 2 _ 2 + Xf 20 d o + X 242 d l 4 _ 4 ] >, (55) 



fjf « — j C Z X S X 22X oorfo2 + Cgl,2— ====(Xu 2 di 1 + Xi3 2 d 3 „!) 



l in in 



where 



2 ^gl,2 [( 2 Xo22XooO + Xf 2 o)<4o + X 2 2oX 2 42d-24] p (56) 



X 022 « e -^ + D-4] CT V4 7 (57) 



X121 
X132 
X242 



_I^(/ + 2 )G-l)e-[ i ( z+1 )- 8 / 3 ] CT2 / 4 , (58) 
-V(' + 3)(l-2)e-i , ( 1 + 1 '- 20 /* 2 / 4 , (59) 



2 



j V(/ + 4)(/ + 3)(Z - 2)(Z - 3)e-W i+1 )- 10 ] ff2 / 4 . (60) 



These expressions for the Xi mn are accurate to 0(a 2 ) at low Z, and have the correct non-perturbative form at high /. 
Since only Xooo and X022 enter at lowest order, the other exponentials may be further safely approximated as ~ Xooo 
since their contributions will be negligible at low I. 
In the limit of zero separation (3 — > we have 

i + (o) = £^(Q*+^) = £^(A B +<5f), (6i) 

i-(0) = 6f(0)=0, (62) 

where Cf and Cf are the lensed E- and -B-mode power spectra respectively. This shows that lensing does not change 
the total polarization power, though it mixes E and B modes as well as different scales. 

2. CMB power spectra 

Once the lensed correlation functions have been computed, transforming to the CMB power spectra is straightfor- 
ward using 

(63) 
(64) 
(65) 
(66) 

For a further discussion of correlation functions and the transform to power spectra see Ref. |l7j . 





cf 


= 2tt 


[ f(Mo(/3)dcos/3, 


Cf 


-cf 


= 2vr 


j £_(/3)4_ 2 (/?)dcos/3, 


cf 


+ cf 


= 2tt 


j £ + (/3)4 2 (/?)dcos/3, 




cf 


= 2tt 


f fx(Mo(/?)dcos/3. 
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FIG. 4: Difference between the lensed and unlensed temperature, cross-correlation and _E-polarization power spectra (top 
three plots), and the lensed B power spectrum (bottom) for a fiducial concordance model. The unlensed model has no tensor 
component (so no B-mode power), and the lensed B power spectrum shown is not highly accurate due to the neglect of 
non-linear evolution in the lensing potential. The magnitude of the lensing effect depends on the fluctuation amplitude in the 
model; here the model has curvature perturbation power A B — 2.5 x 1CP 9 on 0.05 Mpc -1 scales and spectral index n s = 0.99. 



3. Numerical implementation 

The correlation function method is inherently very efficient, only requiring the evaluation of one dimensional sums 
and integrals. For an accurate calculation of Cf it is essential to compute the full range of the correlation function 
because it is sensitive to large and small scales. However, when Cf is not needed the lensing is only a small-scale effect 
and we only need to integrate some of the angular range to compute £(/3) — (and hence the lensing contribution 
Ci — Ci). We find that using /3 max = 7r/16 is sufficient for 0.1% accuracy to I = 2000, providing a significant factor of 
16 gain in speed. Truncating the correlation function does induce ringing on very small scales, so if accuracy is needed 
on much smaller scales the angular range can be increased. For all but very small scales, and the Cf spectrum, we 
can accurately evaluate the sums over I to compute the lensed correlation functions by sampling only every 10th I, 
yielding an additional significant time saving. 

Our code is publicly available as part of CAMB, 3 with execution time being dominated by that required to compute 
the transfer functions for the CMB and the lensing potential. Once these have been computed, the time required to 



3 |http : //camb ■ info] 
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FIG. 5: Comparison of our new result with the O(Cf) harmonic result of Ref. Q (dashed) and the flat-sky non-perturbative 
result of Ref. extended to second order in (7gi,2 (solid). The magnitude of the difference depends on the exact model and 
we have neglected non-linear contributions to the lensing potential. 



compute the unlensed Ci and then lens the result is about a hundred times less (if Cf is not required accurately). 
This means that efficient methods for exploiting 'fast' and 'slow' parameters 0, ^| during Markov Chain Monte 
Carlo parameter estimation can still be applied when lensing is accounted for via the lensed power spectrum. 

Sample numerical results for the lensed CMB power spectra compared to the unlensed spectra are shown in Fig. 0] 



IV. COMPARISON OF METHODS 



We are now in a position to compare our new, accurate full-sky result with previous work. In Fig. [S] we compare 
our result with the full-sky lowest-order perturbative harmonic result of Ref. [correct to 0{Cf)\ for a typical 
concordance model. We also compare to the flat-sky result of Refs. 0,0 which is non-perturbative in a 1 . [We extend 
their results to second order in C g i ; 2 using Eqs. I|17[) and J2J]. In all cases we use an accurate numerical calculation 
of Cf, rather than the Limber approximation, and ignore its non- linear contribution. 

It is clear that the lowest-order perturbative harmonic method of Ref. 0| is not sufficiently accurate for precision 
cosmology, with ~ 1% errors on the temperature and ~ 5% on the E'-mode polarization by I ~ 2000. These errors 
are sufficient to bias parameters even with the planned Planck satellite observations. The perturbative harmonic 



4 |http : //sci . es a . int/planck 
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result is equivalent to expanding the correlation function result self-consistently to first order in Cf . As discussed 
in Sec. 1111 B[ this is inaccurate because l 2 a 2 in the isotropic terms is not very small for large I, so many terms need 
to be retained to get accurate results. It is possible to extend the harmonic result to higher order (2£j, however the 
multi-dimensional integrals required scale exponentially badly with increasing order. Even a self-consistent expansion 
to second order in Cf is not good enough at / > 2000, so at least third order would be required. Furthermore we 
see that for Cj 3 the method is also somewhat inaccurate on large scales: because the B-mode signal comes from a 
wide range of I, and the iJ-mode power peaks on small scales, the non-perturbative effects can be significant on all 
scales. In fact, the large-scale lensed E-mode power also receives most of its contribution from small-scale modes 
since the unlensed polarization power spectrum rises steeply with I on large scales. However, lensing is still only a 
small fractional effect for E'-polarization on large scales and so the perturbative expansion is relatively more accurate 
for E than B. 

The correlation function methods can easily handle the isotropic term non-perturbatively. The accurate flat-sky 
result is much more accurate than the lowest-order harmonic full-sky result, with only ~ 0.3% curvature corrections 
to the temperature. 5 The polarization errors are rather larger, with percent-level difference on Cf ■ Although this is 
smaller than the effect of non-linearities in the lensing potential (see Section^, the latter can be accurately accounted 
for with better modelling (e.g. Ref. |21|) or simulations. While the accurate flat-sky result is probably sufficient to 
Planck sensitivities, curvature effects must be taken account for truly accurate results approaching the cosmic-variance 
limit. Although the curvature is negligible on the scale of the deflection angles, it is not negligible on the scale of 
the lensing potential coherence length. Computing our full-sky accurate result is not much harder or slower than 
computing the flat result, so we recommend our new calculation for future work. 

Note that the absolute precision of the lensed results is limited by the accuracy of the computed lensing potential 
and the unlensed CMB power spectra. In particular, uncertainties in the ionization history may generate errors 
significantly above cosmic variance on the unlensed Ci. We use the RECFAST code of Ref. j2^] that may well be 
inaccurate at above the percent level 6 (2^, [24| . However if the ionization history can be computed reliably to high 
accuracy our new lensing method can then be used to compute the lensed power spectra accurately. 



V. NON-LINEAR EVOLUTION 



The most important assumption we have made so far is that the lensing potential is linear and Gaussian. On small 
scales this will not be quite correct. Although our method does not allow us to account for the non-Gaussianity, we 
can take into account the effect of non-linear evolution on the power spectrum of the lensing potential [and hence 
cr 2 (/3) and C g i,2(/3) Q]- On scales where the non-Gaussianity of the deflection field is small this should be a good 
approximation, assuming we have an accurate way to compute the non-linear power spectrum of the density field. 

We use the HALOFIT code of Ref. [2l| to compute an approximate non-linear, equal-time power spectrum given an 
accurate numerical linear power spectrum at a given redshift. halofit is expected to be accurate at the few percent 
level for standard ACDM models with power law primordial power spectra (but cannot be relied on for other models, 
for example with an evolving dark energy component). We simply scale the potential transfer functions T^{k, rf) of 
Eq. (0 so that the power spectrum of the potential "J has the correct non- linear form at that redshift: 

/-pnon-lincaiYl, 

^^^^fssr ■ (67) 

Since non-linear effects on Cf are only important where the Limber approximation holds, Eq. (|67|l should be very 
accurate. 

The effect of the non-linear evolution on the power spectrum of the lensing potential is shown in Fig. ^ Although 
there is very little effect on scales where the power peaks (I ~ 60), non- linear evolution significantly increases the power 
on small scales. The corresponding changes to the lensed CMB power spectra are shown in Fig. The temperature 
power spectrum Cf is changed by < 0.2% for I ~ 2000, but there are percent level changes on smaller scales. Thus 
inclusion of the non-linear evolution will be important to obtain results accurate at cosmic-variance levels, but is not 
likely to be important at I < 2000 for the near future. The effect on the -B-mode power spectrum is more dramatic, 



5 Due to the opposite sign of curvature and second-order corrections in Cf , the flat-sky correlation result correct to 0{C & \ t 2) is actually 
slightly more accurate than the result correct to 0(C^ 2 ). 

6 http : //cosmocof f ee . inf o/viewtopic ,php?t= 174| 
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FIG. 6: The fractional change in the lensed Ci due to non-linear corrections using halofit |21| for the same model as Fig. 0] 
The lensed Ci are computed using our new accurate method. 



giving a > 6% increase in power on all scales. On scales beyond the peak in the B-mode power (I > 1000) the extra 
non-linear power becomes more important, producing an order unity change in the i?-mode spectrum on small scales. 
On these scales the assumption of Gaussianity is probably not very good, and the accuracy will also be limited by 
the precision of the non-linear power spectrum. For more accurate results, more general models, and on very small 
scales where the non-Gaussianity of the lcnsing potential becomes important, numerical simulations may be required 
(e.g. see Refs. Hill). 

There are, of course, other non-linear effects on the CMB with the same frequency spectrum as the primordial (and 
lensed) temperature anisotropies and polarization. The kinematic Sunyaev-Zel'dovich (SZ) effect is the main such 
effect for the temperature anisotropies, and current uncertainties in the reionization history and morphology make the 
spectrum Cf uncertain at the few percent level at I = 2000 27]. This is a little larger than the error in the first-order 
harmonic lensing result, but this doesn't mean that one should be content with the error in the latter. Precision 
cosmology from the damping tail will require accurate modelling of both lensing and the kinematic SZ effect. Errors 
at the percent level in the lensing power on these scales would seriously limit our ability to constrain reionization 
scenarios with future arcminute-resolution observations. For the polarization spectra, the kinematic SZ effect is much 
less significant [28j. 



VI. CONCLUSIONS 



We have presented a new, fast and accurate method for computing the lensed CMB power spectra using spherical 
correlation functions. Previous perturbative methods were found to be insufficiently inaccurate for precision cosmology, 
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and non-perturbative results in the flat-sky approximation are in error at above the cosmic- variance level. The method 
developed here should enable accurate calculation of the lensing effect to within cosmic-variance limits to I < 2500 
under the assumptions of the Born approximation and Gaussianity of the primordial fields. Non-linear corrections to 
the lensing potential have only a small effect on the lensed temperature power spectrum, but are important on all 
scales for an accurate calculation of the lensed B-mode power spectrum. 
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APPENDIX A: ROTATING SPIN- WEIGHT HARMONICS 

Consider evaluating s Yi m at Dn, where D is the rotation operator corresponding to Euler angles a, (3 and 7. This 
is the same as rigidly rotating the function s Yi m (as a scalar) by Z) _1 and evaluating at n. For spin-0 harmonics we 
know that 

Y lm (Dh) = D l m , m (- % -(3, -a)Y lm ,(h). (Al) 

For spin-s harmonics, we note that 



9/4-1 

s Y lm (n) = (-l)«V_±_zy ms W>,M), (A2) 



where (9, <f>) refer to n, so that 
D 



m 1 in 



(-7,-0,-a) a Y lm ,(n) = (-l) m ' J ^^D l m , m (- % -/?, -a)D l _ m , 0, 0) 



2/ 4- 1 

(-!) ro A/ ft i)dL s (^ e, 0) 



= (-ir^^V ms o^', K ) 

= s Y lm {Dn)e- isK . (A3) 

Here, we have used D(a, (3, y)D((f>, 9, 0) = D{<j>' , 9', k), so that (9', <j>') refer to the image of n under D(a, (3, 7), and k 
is the additional rotation required about Dn to map the polar basis vectors there onto the image of the polar basis 
at n under D(a, (3, 7). Denoting the polar basis (unit) vectors at n by eg and e^, and at n' by e' g and e^, we have 

e' d ±ie^ = e ±lK D(e e ±ie^. (A4) 

This ensures that the 21 + 1 rank-s tensor fields ±Y; m (n) = ± s Yi m (h)(eg + ie l p)<S)- ■ ■ (g) (eg =F* e 0) transform irreducibly 
under rotations as D±Y lm = J2 m ' Dl m >m± Y lm>- 



APPENDIX B: EVALUATION OF A 4 



The integrals 

X imn = J™ ^^Je- a2 ^d l mn (a)da (Bl) 

that are required for the non-perturbative calculation of the lensed power spectra on the spherical sky can easily be 
evaluated as series in a 2 . From the definition of the rotation matrices, we have 

d l mn {a) = {lm\e- iat v\ln), (B2) 



16 



where we adopt the Condon-Shortley phase for the eigenstates \lm) of the L z and L 2 angular momentum operators, 
and we have set H = 1. Expanding the exponential as a series in a, we have 

oo 

X imn = J2 'I U + 1 + i/2){lm\(-iL y y\ln). (B3) 

3=0 J - 

The action of the L y operator on the eigenstates of L z is given by the familiar result 

-iL y \lm) = \ \Jl(l + 1) - m(m - Tj\lm - 1) - \ + 1) - m(m + 1)|Z m + 1), (B4) 

and this can be used recursively to evaluate the matrix element in Eq. (|B3jl . 

For I 3> 1/(7 the series in Eq. (|B3|) is slow to converge and a non-perturbative treatment of the Xi mn is required. 
Using the asymptotic result 

<4» ~ (-l)"- m Jm-n[a + l/2)a], (B5) 

valid for ? — > cx> with S> 1 but a <C 1, and, noting that we only require the case i = m — n, 7 we can use Eq. (6.6314) 
from Ref. [29| to show that 

X mn ~(- { ±±^j e-^ 2 ° 2 '± forz = m-n. (B6) 

In practice, we obtain an excellent approximation to Xi mn , valid for all I, by adjusting the ^-independent term in the 
exponent of the asymptotic result, and the prefactor, so that its series expansion agrees with a direct evaluation of 
Eq. to 0{<j 2 ). 



Note that Xi mn = (—l) rn ^~ n Xi nm so we can always take m> n. 



APPENDIX C: FULL SECOND ORDER RESULT 
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The full result for the lensed correlation functions accurate to second order in C g i and C g \ ; 2 is 8 
£ ~ ^2 47r Cf j (^ooo + 2C g iX oo^ooo + Cgi Woo^ooo + 2A o ) + ;2 ^ooo) ^oo 

+ _j_ (^000 + 2 Cgl^OOo) (Cgl,2^1-1 + Cgl^ll) + ^1,2^220^2-2 

+ Cgl^220^22 ~~ 4CglCgl j2 X 000 X 2 20Ci2 



(CI) 



£+ ~ ^2 ~~ 7 (Cf + Cf) j (^022 + 2C g lX 22^022 + Cgl (J^022-^022 + 2A Q f 2 ) + Cgi 2 A | 2 ) d 22 



Cgl^l32(^132 + 2CgiX{ 32 )(i33 + C g iXi 2 i(Xi 2 i + 2CglX 12 l)rfll 



- 2C g i j2 [X 12 iX 132 + C g \(X 12 iX[ 32 + Xt32^"i2i)] d l 31 + T:C gl X 220 d l m 



+ C gl2 X 22a X 242 d 40 - 2C g \C g i^ 2 (X 220 X 022 d 20 + X 242 X 022 d 42 ) + -C gl X 242 d 44 



i- 4?r ( C F ~ C F)\ (^022 + 2CglX 22^ 22 + Cgl(^022^022 + 2A Q22 2 ) + C^ l2 X' 022 2 ^ d\_ 



(C2) 



V-2Z + 1 x 



+ C g i j2 Xi3 2 (Xi3 2 + 2C g iX 132 )(i3_3 + C g i ! 2-X*i2i (X121 + 2C g \X' l2l )d l _ l 
+ 2C g i (X121X132 + C g i[Xi2iX{ 32 + X 132 X' 12l ]) d l 3 _ 1 + -C gl 2 X 220 d l 00 

+ C^X 220 X 242 d l 40 - 2C g iC g i. 2 (X 220 X' 022 d 20 + X 242 X' 022 d l 4 _ 2 ) + -C^ l 2 X 242 d 4 _ 4 ^ 

(C3) 

^022^000 + Cgl(X022^000 + -^000^022) + Cgl, 2 (-^000^022 + 7^220) 



+ ^ (Jf | 20 + XooO-X"o22 + ^022^000 + 4^022^0 



ooo J 



'20 



2 



[C g i(Xi 2 iX ' 00 + A 121 A 000 ) + Xi 2 iX 000 ] (Cgid^! + C g i )2 dii) 
[C , g i(Xi 32 Xo' 00 + X' 132 X' 000 ) + Xi 32 X 000 ] (C g i i2 d 3 -i + C g irf 31 ) 



C g iCgi i2 (X 000 X 220 (i 00 + ^000^242^40) + -^X 242 X 22 Q(C g ^ 2 d 4 _ 2 + C gl d 42 ) 



(C4) 



As discussed in the main text, the C g i terms may be neglected at the 0(10 4 ) level for realistic lensing deflection 
amplitudes. 
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